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ABSTRACT 


By using the geostrophic approximation in various of the horizontal 
acceleration terms in the Eulerian equations for horizontal flow, wind 
components are computed from which relative vorticity is determined. 
The problem of obtaining a stream function is then reduced to solving 
a Poisson-type equation, thus avoiding the use of the balance equation. 

Two different computation schemes for the x, y components of the 
500-mb wind are investigated. The stream functions obtained from these 
methods are then compared with those obtained through the solution of 
the balance equation. The results of forecasts for 24 and 48 hours 
using the various stream functions with the barotropic model are also pre- 
sented. 

The writer wishes to express her appreciation to Professor George 
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dur ‘Arnason for his many informative discussions and suggestions, and 
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Table of Symbols 
the Coriolis parameter,2 © sin Z, where @ is the 
geographical latitude. 
the mean value of the Coriolis parameter, 2 A. Sin 45°- 
the height of an isobaric surface. 
the upward component of the apparent gravitational 
acceleration. 
the geopotential. $= WK a 
the relative vorticity. °= dv dh 

ax dy 

the absolute vorticity. N= C4£ - 
the stream function for the non-divergent component of 
velocity. 


horizontal Laplacian operator on a constant pressure 


surface. 


iS 


horizontal Jacobian operator, tJ (AB) ss 20 3 = 3h x. 


wt 
Cy 
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Deinwoduction . 

It has been found that the use of the geostrophic approximation in 
the quasi-geostrophic barotropic model for numerical prediction gives 
rise to spurious anticyclogenesis fale . It has further been shown that 
such errors can be eliminated by replacing geostrophic wind and geostro- 
phic vorticity by wind and vorticity derived from a stream function, thus 
retaining only the divergence-free part of the wind 

The present method used by both the Joint Numerical Weather Predic- 
tion Unit (INWPU) in Suitland, Md., and the U. S. Navy Fleet Numerical 
Weather Facility (FNWF), Monterey, Calif., to obtain a stream function 


is through the solution of the balance equation 


a 
NO ae ed en SE 
Tver yvy-evit e Ieee ax | £ (1) 


where ® is the geopotential. This equation, a non-linear partial differen- 
tial equation of the Monge-Ampere type, may be elliptic in one part of the 
map, parabolic in another and hyperbolic ina third area. The problems 
encountered a obtaining significant solutions to this equation have been 
discussed by ‘Arnason [2] : 

The purpose of this research is to examine the possibility of obtaining 
a stream function in a somewhat simpler manner. By using the geostro- 
phic approximation in various of the horizontal acceleration terms in the 
Eulerian equations for horizontal flow, wind components may be computed, 
from which the relative vorticity may be determined. The problem of ob- 


taining the stream function is then reduced to solving the equation 


Seo y (2) 
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a Poisson-type equation readily solved by numerical methods. Thus the 


use of the balance equation may be avoided. 





2. Background. 


Assuming the hydrostatic approximation, and using pressure (p) as 


the vertical coordinate, the equation of motion may be expressed as 


AV = - G¥,Z-£KXV (3) 


In scalar component form, equation (3) becomes, upon expanding 


the acceleration components, 


where w= a ° 

Various approximations have been used to represent the actual wind. 
The simplest is the geostrophic approximation in which the entire accele- 
ration is assumed to be identically zero, yielding 

Va a) f°Y,2 se (5) 

A second approximation widely used is that the flow is gradient. 
With this assumption, tangential accelerations are omitted, but centri- 
petal accelerations are retained. 

The wind law to be used herein primarily consists of omitting the 
local time derivatives. In addition, since the results are to be applied 
to one level only in the present study, the acceleration terms involving 


the vertical velocity advection will be omitted. Thus equations (4) 


reduce to 
YIKaz-=- Mg + fy 
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+tug- fw . (6) 
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By introducing geostrophic wind components in various terms on the 
left-hand side of equation (6) the u, v components of the wind and the 
relative vorticity were obtained. 

Once Aye relative vorticity had been computed, a stream function, 
‘VY , was obtained by solving the Poisson-type equation 

g-L gy (7) 
f: 
using the relaxation technique. The introduction of the coefficient Ee 
is necessary to give the computed W-field units of height. This \Y-field 
as computed was then used to make 24- and 48-hour barotropic prognoses 
which were compared to forecasts of Se the JNWPU and FNWF stream 
functions using the same barotropic model. 

All computations were done on the Control Data Corporation Model 
1604 Computer using the U. S. Navy Fleet Numerical Weather Facility's 
octagonal grid. This grid, centered on the pole and extending to approxi- 
mately 10° N latitude, consists of 1977 grid points with a grid size of 
381 km at 60" latitude. Geostrophic winds and geostrophic vorticity 
were also computed in order that the effects of the approximations intro- 
duced could be compared. Since it is known that the geostrophic approxi- 
mation is poor in low latitudes, and that small errors in the pressure field 
will give relatively large errors in computed winds in this area, the sine 
of the latitude was not allowed to become smaller than sine 15° through- 
out the computations. 

Initial data used were the 500-mb heights for OOZ, 3 January 1958. 
The W-fields for OOZ, 4 and 5 January 1958 were also computed in order 
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to compare the forecast fields with the computed. 





3. Procedures. 
Basically, two computation schemes were investigated to obtain 
approximations of the x, y components of the 500-mb wind. In the first 
method, geostrophic wind components were substituted in the partial de- 

v 
rivatives of equations (6) . Thus the following simultaneous system of 
equations is obtained: 
(1+ d)u 4 (Fe es) v= U. 
Figs) ue (i-3gavey, 


These equations were then solved for u and v, with the following re- 


(8) 


strictions imposed: 

(a) an arbitrary lower limit of .25 was imposed on the terms 
fol aah : us | and (| — F eu ), in order to avoid conditions approaching 
inertial instability; 
and (b) the determinant of the coefficients was not allowed to become 
smaller than .125, to avoid excessive wind speeds. The above method 
will hereinafter be designated as Method I. 

The second method investigated, to be referred to as Method II, 
consisted of replacing all the velocity components in the horizontal ac- 


celeration terms of equations (6) by the geostrophic components. Equa- 


tions (6) then become 


_ Uy dvy _ V5 Ys 
Mealy — “phigh Fo ody 
| yu uD, (9) 
lie wi ose . 


This method, though presumably a cruder approximation to the actual 
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wind, has the decided advantage of avoiding the arbitrary restrictions 
imposed in Method I], as well as a greater simplicity in programming. 

Since the stream field obtained by Method II was very nearly the 
same as that obtained by Method I, one can consider the sequence of 


approximations obtained by the equations 


CEE BW Ch I ae CO 

Valve ee €d aie ae? 
Jf 3x Fo Oy 

(n) in-4) DY (ne) aia) (n-1) (10) 
V am NZ + a CS: aa we Gu 
Vig Ox fF ody 
. (2) (2) 
with n=l-:+°*+ ee  . Method III consisted of computing u andv 


After computing the wind fields for the various methods, and the 
geostrophic wind, =-4 f v,2 XIX, for OOZ, 3 January 1958, the 
difference between computed and geostrophic wind speeds AY | . NAR 
was examined in order to assess the effects of the various approximations 
on the horizontal acceleration. In the regions of sharply defined cyclonic 
curvature, Method II produced the greatest decrease, and, in regions of 
anticyclonic curvature, the least increase in the geostrophic wind. This 
is to be expected as may be seen from gradient-wind considerations. 
When the geostrophic approximation is used in the acceleration terms, 
the latter will tend to be overestimated in regions of cyclonic curvature 
and thus reduce the wind too much below the geostrophic. Similarly, in 
the vicinity of a ridge, the geostrophic approximation in the acceleration 
will tend to underestimate the acceleration and thus the resulting in- 
crease of the wind over the geostrophic will be inadequate. 

The introduction of the higher order approximation by Method III 
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tended to bring the wind speed again toward the geostrophic and, in 
regions of extreme cyclonic curvature, the computed speeds exceeded 
the geostrophic in some cases. This will be further discussed in con- 
nection with vorticity. 

Tables 1 and 2 show representative values of the difference be- 
tween computed and geostrophic wind speeds found in regions of strong 


curvature. 
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Table I; Maximum difference between computed and geostrophic wind 
speeds observed in selected areas of well-defined cyclonic 


curvature. 
Maximum reduction in 
Geostrophic geostrophic wind speed 

Location wind speed Method I Method II Method III 
Eastern 

United States 50 to 58 mps -l17 mps -26 mps -8 mps 
Eastern 

Atlantic 30 to 38 mps -12 mps -19 mps -3 mps 
Eastern 

Europe 30 to 38 mps - 8mps -ll mps -4 mps 
Off East 

Coast of Asia 50 to 62 mps -20 mps -3l mps -9 mps 
Central 

Pacific 50 to 60 mps = Ie Pn ays -28 mps =7 0s 


Table II: Maximum difference between computed and geostrophic wind 
speeds observed in selected areas of well-defined anti- 
cyclonic curvature 





Maximum increase in 


Geostrophic geostrophic wind speed 
Location wind speed Method I Method II Method III 
Off East Coast 
of North America 15 to 20 mps L2smps 7 mps 9 mps 
Siberia 25 to 31 mps 8 mps 6 mps 7 mps 
Central 
Pacific 30 to 39 mps 18 mps ll mps 12 mps 
Western 
Canada 4 to 14 mps 3 mps 2 mps 2mps 
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a 
Next the fields of geostrophic vorticity, = J VY Z , and relative 


Jf 
vortieity, Cz gy ~ se , were computed. As was expected, the de- 


rived values of cyclonic vorticity associated with the migratory lows and 
troughs were smaller than the corresponding geostrophic quantities, with 
the reverse true in the case of anticyclonic vorticity. As shown in Table 
III, Method II produced the greatest decrease in cyclonic vorticity and 
the least increase in anticyclonic vorticity. With the introduction of the 
higher-order approximation, Method III, the values of both the positive 
and negative vorticity centers increased over those obtained from the 
first two methods. Again, in regions of strong cyclonic shear, the maxi- 
mum values exceeded the geostrophic. Further iterations using the tech- 
nique of Method III were made to determine if the introduction of still 
higher-order approximations would indicate convergence. 

The vorticity fields computed from u ey and woe did show some reduc- 


tion in the cyclonic vorticity maxima, but the negative vorticity values 


(4) ) 


kept increasing. With the next iteration, u and v ‘*? | the values ob- 
tained in regions of strong relative vorticity produced overflow in the 
fixed-point computations of this parameter, which had been scaled fora 

-6 =i 
maximum value of 253 x 10 sec. This suggests, perhaps, that the 
sequence of computed winds was not converging toward realistic values, 
at least in one particular region, having given rise to excessive vorti- 
city values. 

After the relative vorticity fields were computed, the Poisson equa- 

4 a . 
tion ee = VY WW was solved for the stream function W using the 
f 
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relaxation technique. This technique as programmed by FNWF scans 
through the field of W laterally using the extrapolated Liebmann method 


of computing a new point value 


Dei 


| ae 

Wie Py ty (11) 
i = So: . 

where Rij is equal to a (Ty, = 5 f) and TF is the relaxation factor, 


in this case 0.35. The 500-mb height field was used as the initial ap- 


proximation. 
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Table II] 


Position of 


associated low 


39N  10E 
58N 22E 
38.5N 33E 
78.5N 88E 
54N 147E 
56N 164W 
84N 135W 
55N 71W 
46N 23W 


67 


80 


80 


38 


73 


09 


or 


84 


100 


Maximum value of center of 
cyclonic relative vorticity 


Method I 


52 


68 


63 


34 


61 


Sell 


ell 


79 
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Method JI 


47 


63 


08 


34 


o7 


50 


90 


71 


og 


Method JJI 


60 


70 


73 


36 


65 


yl 


o2 


74 
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4. Results of Stream-Function Computations. 

A comparison of the computed stream-function fields with those ob- 
tained by the FNWF solution of the balance equation reveals that, in all 
instances, the gradient of Y between the pole and the grid boundaries is 
slightly greater for the fields computed by the methods of this study than 
for those computed through the balance equation, the difference varying 
from 200 to 400 feet in the cases studied Through the relationship 
WV =-YW xix , it is obvious, then, that on the average the non- 
divergent part of the wind, computed by the methods described herein, 
slightly exceeds that obtained through the balance equation. Since, 
however, the difference in the gradient is almost uniform from the bound- 
aries to the pole, the difference in wind speed at any point would be only 
of the order of centimeters per second. 

Considering only the stream-function fields computed by the methods 
of this study, the gradient of W between pole and grid boundaries was 
least for Method II, corresponding to the lowest wind speeds. 

A comparison of the computed stream functions with the FNWF 
stream functions and the height fields themselves reveals that the place- 
ment of the major systems is the same in all cases. However, an exam- 
ination of the stream fields obtained by Method I reveals that the decrease 
in the amplitude of the troughs is slightly greater than that observed be- 
tween the height and the FNWF stream-function fields. This decrease 
in amplitude of the major troughs is extremely pronounced in the 
\/-fields computed by Method II, presumably due to the underestimation 
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of the wind by this method as previously mentioned. 

With the introduction of the higher order approximation, Method III, 
this pronounced weakening of the troughs is not observed, the amplitude 
of the features being approximately the same as those of the FNWF 


\ -fields ; 
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5. Results of Stream Barotropic Prognoses. 

Using the stream functions obtained by the various methods, 24- 
and 48-hour prognoses were made using the barotropic forecasting equa- 
tion 


23W sit De 
vst +o} (yn) : = 6 


tle 


(12) 


where aif is the Jacobian operator (3.] . The term — An wv is added to 
y gt 

the non-divergent barotropic model to stabilize the ultra-long atmos- 

pheric waves. In these computations _. was assigned the value of 4, 

in accordance with the results of Cressman [4]. 

In order to have some basis for comparison, 24- and 48-hour prog- 
noses were also made using the same model with the stream functions 
obtained through both the FNWF and the JNWPU solutions of the balance 
equation. For verification purposes, difference fields were formed and 


the "pillow" and root-mean-square error computed over the entire grid 


using the formulae 


197T 


= (A-B),, 
Pillow = ™—-—— (13) 
Cie Le ae 
| (A-8)- Prcow | 
RMSE =7- 1) el (14) 


where A is the computed stream function value and B is the forecast 
value. 
As can be seen from Table IV, no significant difference in either 


the “pillow” or root-mean-square error is apparent for the case studied. 
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However, it may be noted that the root-mean-square error for both the 
24- and 48-hour forecasts is least for Method II and largest for Method 
II. 

An examination of the error fields themselves showed a marked simi- 
larity in the regions of errors greater than 200 feet. The only major dif- 
ferences that could be noted were in an area north of the Hawaiian 
Islands, where Method III had a much larger error than other methods, 
and in an area south of Greenland, where both Methods I and II produced 
errors approximately 200 feet less than the other methods. 

Since the barotropic model used for forecasting does not include 
provisions for development, only the movement of the major systems was 
used to compare the various methods. Considering only this aspect, no 
great differences appeared with regard to most of the major systems. The 
only notable exceptions were the troughs in the Eastern Pacific and the 
Eastern Atlantic where movement was overestimated less using the 
FNWF stream function than by the other methods. In addition, there was 
a slight indication that, at least in this case, the systems tended to 


move slightly faster when Method I was used. 
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Table IV 


Method used to obtain 24-hour baro- 48-hour baro- 
initial stream function tropic forecast tropic forecast 
Pillow RMSE Pillow RMSE 


FNWF solution of 


balance equation + 8 ft. 174 it. + 6 ft. Vd ie A 


JNWPU solution of 


balance equation -36 ft. 172 ft. -70 ft. ZO bss 
Method ] -14 ft. VPoeorrt. -56 ft. 202. 
Method II +9 ft. 164 ft. +24 ft. 254 fit 
Method III -17 ft. 182 ft. -62 ft. oi 1 fee 
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6. Conclusions. 

Since the methods described herein have been applied to only one 
synoptic situation in this study, only tentative conclusions can be made 
at this point as to the relative merits of the various ways of obtaining a 
stream function for use in numerical prognostic models. All methods 
have shown comparable results and, until many different weather pat- 
terns are tested, it is difficult to separate the errors in the forecasts 
into those arising’ from inadequacies in the barotropic model and those 
due to the assumptions made in the various methods. 

It has been seen, however, that the stream-function fields obtain- 
ed from Method II had the least computed root-mean-square error of any 
of the methods tried anda "pillow" smaller than those computed using 
the other two methods of this study or the JNWPU solution to the balance 
equation for both the 24-and the 48-hour prognoses. However, a greater 
decrease in amplitude in the transition from the height to the stream- 
function field has also been observed when this method was used and 
must also be considered in evaluating this method. Although this method 
presumably underestimates the total wind speed, in the case studied the 
24-hour prognostic movement of the major systems did not differ appre- 
Cciably from that of the other methods used. Since no arbitrary restric- 
tions need be imposed when this method is used, and since it is easier 
to program, further study is certainly warranted. 

Since the introduction of the higher-order approximation into the 
same basic computation scheme (Method III) produced spurious results 
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in both the wind and the relative vorticity fields, the use of this pro- 
cedure appears questionable. However, the introduction of certain re- 
strictions on the values of a and a similar to those imposed 
in Method I, may possibly eliminate the difficulties encountered. 

The procedure described under Method I apparently has some merits, 
although further study should be made of the arbitrary restrictions placed 


on this method for the computations of this study. It would also be most 


interesting to try further iterations of this method using equations of the 


form . , joa : yum 
use juve “Fe ) “2 Ug 
au”) Nn (15) 
es?) urs ( — “)v"aNg 


to determine what the effect on the geostrophic deviation and relative 
vorticity would be in comparison to the effects observed upon introduc- 
tion of the higher approximations in Method III. Since extensive repro- 
gramming would be involved in doing this, such iterations were not 


attempted in this study. 
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Appendix I] 


The following charts were computed during the course of this in- 


vestigation and are available upon request from the U.S. Naval Post- 


graduate 
(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 
(10) 
(11) 
(12) 
(13) 


(14) 


(15) _ 


(16) 
(17) 
(18) 
(19) 


(20) 


School, Monterey, California. 

500-mb heights, OOZ, 3 Jan. 1958 

500-mb heights, OOZ, 4 Jan. 1958 

900-mb heights, OOZ, 5 Jan. 1958 

Geostrophic Wind Speed, OOZ, 3Jan. 1958 
Method I, 500-mb Wind Speed, OOZ, 3 Jan. 1958 
Method II, 500-mb Wind Speed, OOZ, 3 Jan. 1958 
Method III, 500-mb Wind Speed, OOZ, 3 Jan. 1958 
Geostrophic Relative Vorticity, OOZ, 3 Jan. 1958 
Method I, Relative Vorticity, OOZ, 3 Jan. 1958 
Method I], Relative Vorticity, OOZ, 3Jan. 1958 
Method III, Relative Vorticity, OOZ, 3 Jan. 1958 
FNWF Stream Function, OOZ, 3 Jan. 1958 

Method I, Stream Function, OOZ., 3 Jan. 1958 
Method II, Stream Function, OOZ, 3 Jan. 1958 
Method III, Stream Function, OOZ, 3 Jan. 1958 
FPNWF Stream Function, OOZ, 4 Jan. 1958 

Method I, Stream Function, OOZ, 4 Jan. 1958 
Method II, Stream Function, OOZ, 4 Jan. 1958 
Method III, Stream Function, OOZ, 4 Jan. 1958 


FNWF Stream Function, OOZ, 5 Jan. 1958 
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(21) 
(22) 
(2 3) 


(2 4) 


(25) 


(26) 


(27) 


(28) 


(29) 


(30) 


(31) 


(32) 


Method I, Stream Function, OOZ, 5 Jan. 1958 
Method II, Stream Function, OOZ, 5 Jan. 1958 
Method III, Stream Function, OOZ, 5 Jan. 1958 


FNWF 24-hr. Stream-Function Prognosis, Verifying Time, 
OOZ, 4 Jan. 1958 


Method I, 24-hr. Stream-Function Prognosis, Verifying 
Time, OOZ, 4 Jan. 1958 


Method II, 24-hr. Stream-Function Prognosis, Verifying 
Time, OOZ, 4 Jan 1958 


Method III, 24-hr. Stream-Function Prognosis, Verifying 
Time, OOZ, 4 Jan. 1958 


Difference between FNWF Computed W‘\-field and 
Prognosis, OOZ, 4 Jan. 1958. 


Difference between JNWPU Computed \/-field and 
Prognosis, OOZ, 4 Jan. 1958 


Difference between Method 1 Computed \) -field and 
Prognosis, OOZ, 4 Jan. 1958 


Difference between Method II Computed \)-field and 
Prognosis, OOZ, 4 Jan. 1958 


Difference between Method III Computed \-field and 
Prognosis, OOZ, 4 Jan. 1958 
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